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Intracellular transport of large cargoes, such as organelles, vesicles or large proteins, is a complex 
dynamical process that involves the interplay of ATP-consuming molecular motors, cytoskeleton 
filaments and the viscoelastic cytoplasm. The displacements of particles or probes in the cell cy- 
toplasm as a function of time are characterized by different (anomalous) diffusion regimes. We 
investigate here the motion of pigment organelles (melanosomes) driven by myosin-V motors in 
Xenopus laevis melanocytes using a high spatio-temporal resolution tracking technique. By analyz- 
ing the mean square displacement (MSD) of the obtained trajectories as a function of the time lag, 
we show that the melanosomes display a transition between subdiffusive to superdiffusive behavior. 
A stochastic theoretical model is introduced to generalize the interpretation of our data. Starting 
from a generalized Langevin equation that explicitly considers the collective action of the molec- 
ular motors we derive an analytical expression for the MSD as a function of the time lag, which 
also takes into account the experimental noise. By fitting our model to the experimental data we 
were able to discriminate the exponents that characterize the passive and active contributions to 
melanosome dynamics. The model also estimates the "global" motor forces correctly. In this sense, 
our model gives a quantitative description of active transport in living cells with a reduced number 
of parameters. 

PACS numbers: 87.16.-b, 87.16. Uv, 87.10.Mn, 87.80.Nj, 87.16. Nn 



I. INTRODUCTION 

Molecular motors are responsible for the active trans- 
port of organelles and other cargoes along cytoskele- 
ton tracks to their correct destination in the cytoplasm. 
Three different classes of molecular motors are involved 
in this task: dynein and kinesin motors, which trans- 
port cargoes toward the minus and plus ends of micro- 
tubules, respectively, and myosin motors, responsible for 
the transport along actin filaments toward the barbed 
end (reviewed in While properties of microtubule- 
dependent transport, considered to be responsible for 
long-distance transport, have been extensively studied in 
vivo 0, H, 0, IE] J we still do not know important aspects of 
the properties and regulation of actin-dependent trans- 
port which are believed to support local, short-distance 
movement of cargoes in living cells Q . 

The cytoskeleton (CSK) is a crowded network of semi- 
flexible linear protein polymers with a complex dynamics, 
that can exert forces and affect the rheology of cells 0, Q . 
Particularly, F-actin network has no global directional- 
ity and consists on intercepting randomly distributed fil- 
aments 9]. The cytoskeleton, organelles, proteins and 
motors are part of the highly crowded cytoplasm where 
the transport of cargoes takes place. 
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Active (AMR) and passive (PMR) microrheolgy are 
two distinct approaches to study the mechanical prop- 
erties of the intracellular environment [lOj. While in 
AMR micrometer-sized embedded probe particles are 
manipulated by external fields and their displacement is 
measured, the spontaneous displacement fluctuations of 
probe particles is analyzed in PMR. 

Recent studies using AMR and PMR have shown that 
the action of ATP-consuming molecular motors drives 
the system out of equilibrium. The violation of the fluc- 
tuation dissipation theorem (FDT) has been observed 
in an in vitro system consisting of a cross-linked actin 
network with embedded myosin motors and also for 
beads attached to the cell membrane of human airway 
smooth muscle (HASM) cells Likewise, direct evi- 
dence of the deviation from equilibrium inside the cell 
has been recently obtained [ll|, [l^. In Ref.p^ both 
forced and spontaneous motions of magnetic microbeads 
engulfed by Dyctyostelium cells are analyzed to derive the 
power spectrum of forces acting on intracellular phago- 
somes. However, the theoretical model proposed to de- 
scribe their data did not distinguish between active and 
passive contribution to the transport. 

In most of PMR experiments it has been observed that 
the mean square displacement (MSD) of the probes dis- 
plays a crossover between a subdiffusive and a superdif- 
fusive regimes ^8, 13, 14, 15, 16, 17, 18]. The subdiffusive 
behavior is characterized by an exponent ranging from 
0.2 [Hi to < 1 0, mi while the superdiffusive behav- 



ior presents exponents close to 1.5 @, Typically, the 
transition time between these two regimes is on the order 
of 1 s. 

While some authors attribute the subdiffusive behavior 
to elastic trapping |2fli] , obstruction |[22|] , crowding [23!] or 
stalling P , others propose that apparent subdiffusion can 
arise from noise inherent to single particle tracking (SPT) 
experiments due to slight errors on the determination of 
the actual particle position Although there is not 

a general consensus of the causes of subdiffusion yet, it 
is well accepted that superdiffusion has its origin in the 
collective behavior of molecular motors 0, [ill, . 

Properties of motor proteins have been well character- 
ized in vitro but not in the viscous environment of a cell. 
Recently, pioneering works developed new techniques 
that allow to measure the forces exerted by molecular 
motors directly [1^, [2^ . In Ref. [IHl motion enhanced 
differential interference contrast (MEDIC) movies of liv- 
ing NT2 (neuron-committed teratocarcinoma) cells at 
37C, was used to determine the force-velocity curves in 
vivo. In Ref. ^2§\ the force that kinesin-1 exerts on lipid 
droplets in fly embryos, was measured using a novel de- 
velopment of optical traps that can be used to apply 
precise forces to moving droplets in vivo. Unexpected 
differences between motor regulation in vivo and in vitro 
were found, revealing that further investigation should 
be done in the field. 

Melanophore cells are an exceptionally convenient 
model system to study intracellular transport driven 
by molecular motors [27|, and thus to investigate some 
properties of out-of-equilibrium systems in living cells. 
Melanocytes have pigment organelles called melanosomes 
which are filled with the black pigment melanin. Then, 
they can be easily imaged using bright-field transmis- 
sion light microscopy (Fig[T|) and tracked with millisec- 
ond temporal resolution and nanometer precision _3]. Us- 
ing a single particle tracking (SPT) technique 0, ll^l, we 
follow the motion of myosin-V driven melanosomes along 
actin filaments and compute the MSD from the analysis 
of their trajectories. 

Melanosomes are more or less spherical and stiff par- 
ticles with sizes around 500 nm. Pigment organelles can 
be distributed in the cells in two configurations: either 
aggregated in the perinuclear region or homogeneously 
dispersed in the cytoplasm. The transport of pigment 
organelles during aggregation and dispersion is regulated 
by signaling mechanisms initiated by the binding of spe- 
cific hormones to cell surface receptors [1^. In a recent 
paper, we analyzed the MSD dependence with the time 
lag using an empirical model [l3[. With the aid of nu- 
merical simulations, we concluded that the main differ- 
ence between aggregation and dispersion condition was 
the average time the melanosome spends diffusing during 
the intervals between successive processive motions. As 
a consequence, active part of the transport would not be 
influenced by stimulation conditions. 

In this work we explore an alternative and complemen- 
tary approach to obtain quantitative information about 
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FIG. 1: A Xenopus laevis melanophore image obtained us- 
ing bright-field transmission light microscopy. Melanosomes 
can be clearly identify and followed with high spatio-temporal 
resolution using a SPT technique. 

the intracellular transport, that can be generalized to 
other biological data. Starting from a stochastic model 
that considers the viscoelasticity of the intracellular envi- 
ronment and the action of molecular motors explicitly, we 
derive an analytical expression for the MSD as a function 
of the time lag, which also takes into account the experi- 
mental noise. By fitting the logarithmic derivative of the 
MSD-versus-lag time to the experimental data we were 
able to discriminate the passive and active contributions 
to the melanosome dynamics. Experiments using cells 
expressing a dominant negative construct of myosin-V 
where melanosomes are not being actively transport by 
myosin-V, confirmed the robustness of the model. Fur- 
thermore, our approach enables us to quantitatively as- 
sess the motor forces correlation function, to estimate the 
magnitude of motor forces and to determine an effective 
diffusion coefficient by only using data coming from a 
PMR experiment. 

II. MATERIALS AND METHODS 
A. Melanophore cell culture and transfection 

Immortalized Xenopus laevis melanophores were cul- 
tured as described in Ref. [3Q[ . In order to track 
the movement of individual organelles, the number of 
melanosomes in the cell was reduced by treatment with 
phenylthiourea 

To study transport along actin filaments, the cells were 
incubated at 0°C for 30 min with 10 /iM nocodazole to 
depolymerize microtubules [111 . 

Melanophores were stimulated for aggregation or dis- 
persion with 10 nM melatonin or 100 nM MSH, respec- 
tively. The samples were observed between 5 and 15 min 
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after stimulation. All the measurements were performed 
at 21°C. 

Cells were transfected using the FuGENE 6 transfec- 
tion reagent (Boehringer Mannheim Corp.) following 
the vendor's protocols with a plasmid encoding a green 
fluorescent protein-tagged myosin-V tail. Expression of 
this plasmid results in a dominant-negative inhibition of 
myosin-V driven melanosome transport '31]. The plas- 
mid was a kind gift of Dr. Vladimir Gelfand (Northwest- 
ern University, Chicago, IL) 



B. Samples preparation for imaging 

For microscopy measurements, cells were grown for 2 
days on 25- mm round coverslips placed into 35-mm plates 
in 2.5 ml of the medium. Before observation, the cover- 
slips were washed in serum-free 70% L-15 medium and 
mounted in a custom-made chamber specially designed 
for the microscope. 



C. Tracking experiments 

Single particle tracking experiments of melanosomes 
moving along actin filaments in wild type cells were car- 
ried out in a Zeiss IM 35 microscope adapted for SPT us- 
ing a 63x oil-immersion objective (numerical aperture = 
1.25) under illumination with a tungsten-halogen lamp. 
A CCD camera (PixelVision, Oregon, US) was attached 
to the video port of the microscope for imaging the cells. 
Movies were registered at a speed of 14 frames/s. 

Tracking experiments of melanosomes in transfected 
cells were carried out in an Olympus 1X70 microscope 
using a 60x water-immersion objective (numerical aper- 
ture — 1.2). A cMOS camera (Pixelink, Ottawa, Ontario, 
Canada) was attached to the video port of the microscope 
for imaging the cells at a speed of 50 frames/s. 

Trajectories of melanosomes were recovered with 2 
nm precision and 10 ms temporal resolution from the 
movies registered as described above using the pattern- 
recognition algorithm described in [3^ . This algorithm 
is included in the program Globals for Images devel- 
oped at the Laboratory for Fluorescence Dynamics (UCI, 
Irvine, CA). The program, which also contains some of 
the tools used for trajectory analysis, can be downloaded 
from the Laboratory for Fluorescence Dynamics website 
(www.lfd.uci.edu) . 



III. EXPERIMENTS 

Melanophores were treated with nocodazole as de- 
scribed in Materials and Methods in order to depoly- 
merize microtubules. After this treatment, aggregation 
and dispersion of melanosomes were induced by addition 
of melatonin and MSH, respectively. Movies of regions 




FIG. 2: (color online) A representative trajectory of a 
melanosome followed during 70 s with 0.07 s resolution. The 
lighter colored segments represent the periods of rectilinear 
motion where the speed of the motor-driven motion was com- 
puted. The cell was stimulated with 100 nM of MSH. Scale 
bar = 0.5/im 

of the cells were recorded from which a total of 134 tra- 
jectories of melanosomes moving along actin filaments 
were obtained in aggregation and dispersion by using the 
pattern-recognition algorithm (3^ . 

Figure [2] shows a representative trajectory of a myosin- 
V driven melanosome during dispersion. The tracks are 
noisy; however, periods of rectilinear motion can be dis- 
tinguished from periods of random or diffusive behavior. 
Analysis of trajectories during the rectilinear segments 
gives a rough estimation of the speeds distribution that 
can be associated with the motor driven motion. We ob- 
tained a wide distribution for the speed v that can be 
well described with a lognormal distribution with mean 
value 70 ± 17 nm/s and standard deviation of log v equal 
to 0.6 ± 0.2 [l^. These values are in agreement with 
the ones mentioned in the literature for myosin-V in vivo 
P|. The noise of the experimental trajectories, quanti- 
fied using a smoothing procedure was between 1 x 10~^ 
and 18 x 10~^ /im^, which is much larger than the noise 
expected from the tracking method (^ 0.2 x 10~^ fj,m^) 

The mean square displacement (MSD) for every tra- 
jectory is calculated as follows, 

MSD{t) = (^{x{t + r) - x{t)f + {y{t + r) - y{t)f) (1) 

where x and y are the coordinates of the particle, r is a 
lag time and the brackets represents the time average . 

Figure[3]shows the average value of MSD obtained dur- 
ing dispersion and aggregation. The average distance 
traveled by the organelles during dispersion is signifi- 
cantly higher than during aggregation as was previously 
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tively transported by myosin-V, we transfected Xenopus 
nielanophores cells with a plasmid encoding a green fluo- 
rescent protein-tagged myosin-V tail. Since this mutant 
of myosin-V cannot attach to actin filaments, expres- 
sion of this plasmid results in a dominant-negative in- 
hibition of myosin-V driven melanosome transport [3l| . 
The MSD(t) was calculated for every trajectory and the 
average behavior is shown in Fig. [31 We also analyzed 
the value of MSD obtained by extrapolation at r = 
s for the dominant-negative cells and verified that the 
MSDo values obtained followed a Gaussian distribution 
with 0.56 X 10"'* iim? mean and cr = 0.15 x 10"'' fim?. 
These values are in agreement with the ones expected 
according to the error of the tracking method [l^, H^l . 

IV. THEORETICAL MODEL 



FIG. 3: (color online) Double-logarithmic plot of the mean 
square displacement (MSD) as a function of the time lag. Ex- 
perimental data points for wild-type cells are symbolized by 
squares (aggregation) and circles (dispersion), while the data 
corresponding to mutant cells are symbolized by rhomboids. 
The symbols represent an average over all the trajectories. 
Inset: MSD vs time lag plot for all trajectories in wild-type 
cells. 

observed [1, HI. For times larger than 1 second, the be- 
havior presents anomalous diffusion with an exponent 
around 1.35. 

The observation of a residual value of the MSD at 
T (MSDo) suggests that transported melanosomes 
move in a constrained region with a dynamics faster than 
the temporal resolution of the tracking method. The 
presence of this jittering also results in an a ppa rent sub- 
diffusive behavior for short times scales [2^ l33j. 

We compute MSDo by linear extrapolation of the first 
four data points of each experimental MSD vs r plot. The 
obtained MSDo values followed an exponential distribu- 
tion, with characteristic value equal to 13±6 x lO"** jirri^ 
for aggregation and 2.6 ± 0.1 x 10~^ fjim? for dispersion. 

It is important to mention at this point that the MSD 
vs r plots for the different trajectories present a wide 
distribution, principally related to the spreading of the 
MSDo values (see the inset in Fig 13). This observation 
was also found in other systems p!3. f33|. In this sense, 
the ensamble average plot shown in Fig. [3] represents 
a poor estimator of the behavior of the driven particle. 
Instead, the logarithmic derivative of the MSD-versus-lag 
time /3(r), defined as [l3| 
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redounds on a more accurate magnitude to describe the 
dynamics of the transport. Figure|4]shows /3 as a function 
of the time lag for the different experimental conditions. 

In order to characterize the dynamical proper- 
ties of melanosomes when they are not being ac- 
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diffusive behavior of a melanosome 
in the intracellular medium and si- 
by molecular motors. For this pur- 
resulting dynamics by means of the 
equation (GLE) 



mX{t)+ dt' -f{t-t')X{t') ^ F{t) 



(3) 



where "f{t) is the dissipative memory kernel that charac- 
terizes the viscoelastic properties of the cytoplasm. The 
random force F{t) is assumed to be the sum of two un- 
correlated contributions 



F{t) = at) + x{t) 



(4) 



being ^(i) the internal noise due to thermal activity, 
while x{t) is an external random force representing the 
action of the molecular motors. 

The thermal noise ^(t) is a zero-centered and station- 
ary random force with correlation function 



{amt')) = c{\t-t'\) 



(5) 



and is related to the memory kernel 7(t) via the 
fluctuation-dissipation theorem [36| 



C{t) = kBTj{t), 



(6) 



where T is the absolute temperature, and fc^ is the Boltz- 
mann constant. This internal noise is the responsible for 
the passive transportation. To reproduce the subdiffu- 
sive behavior due to the viscoelastic properties of the 
cytoplasm, the noise autocorrelation function ([S]) can be 
modeled as a power- law 37 1 



C{t) 



Co 



r(i-A) Vto 



t 



(7) 



where < A < 1, Cq is a proportionality coefRcient, tq is 
set to 1 s., and T{z) is the Gamma function. 
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On the other hand, the random force x{t) is not re- 
lated to the memory kernel j{t) because it stems from 
the activity of the molecular motors. The irreversible 
conversion of chemical energy from ATP hydrolysis into 
the particle motion via the activity of myosin-V motors 
drive the system to an out-of-equilibrium situation. As 
a consequence, the FDT is no longer valid. 

Assuming that the actin network has no global direc- 
tionality, the random force x(0 exerted by the action 
of the myosin-V motor is chosen as a zero-centered one. 
Moreover, in recent experiments it was observed that 
the forces power spectrum exhibits a power law behavior 
BOM 111- Then, we assume that the autocorrelation 
function 



(x(t)x(0>=A(|i-t'|) 



can be written as 13£ 



Ao 



r(l-a) V^o 



(8) 



(9) 



where < a < 1, Ag is a proportionality coefficient and 
To is set to 1 s. 

Considering that the average melanosome diameter is 
500 nm [s^l, its density is ~ 1.2 g/ml [H and 7 ~ lO'^ 
Ns/m (see Refs. in [33]), the characteristic time would 
be ~ 10~^^ s. Then, from an experimental point of view, 
the inertial effects can be neglected. In this situation, the 
temporal evolution of the displacement can be obtained 
from ([3]), which is given by 

Xit) = X0+ f dt'G(t~t'mt')+x{t')), (10) 



where xq — X{t — Q) is the deterministic initial position 
of the particle. The relaxation function G{t) is the inverse 
form of the Laplace transform 



G{s) = 



1 



si{s) ' 



(11) 



where 7(5) is the Laplace transform of the dissipative 
memory kernel. 

To obtain an analytical expression of the MSD it is 
necessary to calculate the two-time position correlation. 
Starting from Eq. (fTO|) and using relation (l6|), it can be 
written as 

{X{t + T)X{t)) = xl + keT {I{t) + I{t + r) - I{t)) 

+ / dh {G{h)H{h + r) + G{h + T)H{h)) , 
Jo 

(12) 

where the involved relaxation functions are given by 



m 



dt'G{t') 



H{t) = f dt'G{t')A{t-t'). 
Jo 



(13) 
(14) 



I{t) only depends on the internal noise, while H{t) is 
related to both internal and external contributions. 

Using the autocorrelation functions ([7]) and (HI), the 
involved relaxation functions can be written as 



m - 



G{t) 



knT 1 



Co r(A+l) \toJ ' 
ksT 1 / t ^ 



roCor(A) V'^o 



H{t) = ekeT 



r(A-a + l) V^o 



where 



Ao 
Co 



(15) 
(16) 
, (17) 

(18) 



is a dimensionless parameter that measures the relative 
intensity among the motors forces and the thermal forces. 

For 2A — a > the integral term of (|12p can be explic- 
itly evaluated. In this case, it can be demonstrated that 
the MSD long time limit 



MSD{t) = lim {{X{t + r) - X{t)Y) 

t^OO 



(19) 



can be written as [3f 



MSD{t) 



^kgT 
70 



r(A + i)Vo 



(-)^-H£ifA,a(-)'^-" 



To 



where 



i^A,a = r(a - 2A) 



sin(7r(A — a)) — sin(7rA) 



(20) 



(21) 



is a positive constant and 70 = Co/kT. The two- 
dimensional situation is considered by the factor 4. 

It is worth mentioning that, although < a < 1, when 
1 < 2A - a < 2 the MSD (gOl) exhibits a crossover from a 
subdiffusive regime with an exponent A, to a superdiffu- 
sive regime with an exponent 2A— a [35| . The first term of 
(PO)) represents the passive transport of the melanosome 
in the viscoelastic medium, while the second one corre- 
sponds to the contribution of the random force x{t) ^^^d 
is originated in the activity of the myosin-V motors. 

Finally, to compare the analytical expression ([^0]) with 
the experimental data obtained using SPT techniques, 
it is necessary to take into account the error on the 
particle position determination. As is established in 
Refs.[23.l40|. the effect of an uncorrelated noise of vari- 
ance 77^ -generated by measurement errors in particle lo- 
cation or by biological activity- is to add a constant to 
the mean square displacement. Then, the MSD{t) given 
in (I^D]) rewrites as 



MSD{t) MSD{t) + {2rif 



(22) 
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As it was previously discussed, the corresponding local 
slope of the MSD-versus-lag time ([2]) is a better function 
to characterize the dynamics. From ((20|) and ([22|l we get 



r(A+i) V TO ' 



r{\+i)yTo' -1-0 



where 



<5 = 70 



(2ry)2 



In particular, setting e = in Eq. (|23p one gets 



^ TO ' 



(23) 



(24) 



(25) 



which corresponds to the absence of active transport case. 



V. RESULTS AND DISCUSSION 

In this paper we study the motion of 500 nm 
melanosomes in Xenopus laevis melanophores treated 
with nocodazole (a microtubule depolymerizer), where 
the active transport is powered by the actin- 
dependent motor myosin-V. Aggregation and dispersion 
of melanosomes were induced by the addition of mela- 
tonin and MSH, respectively. 

The MSD of melanosomes in both stimulation condi- 
tions displayed a transition from a subdiffusivc to a su- 
perdiffusive regime and the corresponding local slope of 
the MSD-versus-lag time was well described by Eq. ([23]) . 
as shown in FigHl 

Four parameters characterize the behavior: A, a, e and 
6, where A and a are the power law exponents of the 
internal and external noise correlation functions, e is a 
parameter that measures the relative intensity between 
random forces and S is associated with the residual value 
of the MSD as r ^ 0. 





A 


a 


e 


S 


aggregation 


0.96 ± 0.04 


0.58 ± 0.08 


202 ± 38 


41 ± 8 


dispersion 


0.98 ± 0.03 


0.58 ±0.10 


83 ±13 


4.0 ±0.6 


mutant m-V 


0.94 ± 0.04 






0.061 ±0.014 



TABLE I: Dimensionless parameters of the model A, a, e and 
S were obtained from the fit of Eqs. 1231 and to the exper- 
imental data shown in Fig. ID The values are given as mean 
± standard deviation of the ensemble. 



The parameters A , a , e and S obtained from the 
fitting shown in Fig [4] are displayed in Table HI The 
power-law exponents A and a for both dispersing and 
aggregating melanosomes were similar (A 0.97 and 
a ~ 0.58). These values lead to an asymptotic exponent 




FIG. 4: The logarithmic derivative of the mean square dis- 
placement as a function of the time lag (Eq. [2]|. The exper- 
imental data for wild-type cells are symbolized with squares 
(aggregation) and circles (dispersion). The solid lines repre- 
sent the fit to P of the expression given by Eq. ((23}. Inset: 
The same plot for mutant cells. The solid line represents the 
fit of Eq. (|25[) . The symbols represent ensemble averages. 



2A— a ^ 1.36, which is a signature of superdiffusion. Sim- 
ilar values were obtained in other systems [1, [13, [13, [l^ ■ 
On the other hand, the obtained value a ~ 0.58 implies 
a power spectrum of fluctuating motor forces that scales 
as oj~°'^'^, corresponding to a smoothing of discontinu- 
ities in instantaneous force pulses, as was stated in Refs. 

S[il. 

It is worth pointing out that with our model we are able 
to discriminate the contributions of internal and exter- 
nal forces in a straightforward fashion. In a recent paper 
[T^ l , Wilhelm described the dynamics of forced motions 
of magnetic microbeads in Dictyostelium cells in terms 
of a generalized Langevin equation. However, differently 
of what we have done in this paper, that model consid- 
ers a single force term which included contributions of 
both thermal Brownian forces and driving forces gener- 
ated by molecular motors. The force correlation function 
followed a power-law behavior with exponent 0.8, which 
is between the values of A and a obtained by us. 

Melanosomes in cells expressing a dominant-negative 
myosin-V construct, and thus not being actively trans- 
ported by myosin-V, show spontaneous motion and the 
MSD displays a subdiffusivc behavior. In the absence 
of motors, Eq. (j23p rewrites in terms of A and 6 (Eq. 
([25)1 ). Interestingly, we obtained a value for A ^ 0.94, 
very similar to the one obtained in the presence of ac- 
tive motors. This reinforces our hypothesis that A is the 
parameter that characterizes the subdiffusivc motion of 
the melanosomes in the cytoplasm. Therefore, our model 
allows us to discriminate between passive and active con- 
tributions to the motion. In both cases (i.e. in the pres- 
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ence and in the absence of active motors), A is close to 1. 
Similar values were obtained for the subdiffusive anoma- 
lous exponent of tracer proteins in highly concentrated 
random-coil polymer and globular protein solutions that 
mimic the crowded conditions encountered in cellular en- 
vironments [23l |. and for the anomalous diffusion of dcx- 
tran polymers inside HeLa cells 21]. 

Typically, the anomalous diffusion regime is character- 
ized by fitting an empirical power-law relationship to the 
MSD dependence with t of the form [H, [H, [4l| : 



Using definitions ((T8)l and (124]) . Aq rewrites as 



MSD{t) = c + D*{ — y, 

To 



(26) 



being D* an effective diffusion coefficient and 7 an expo- 
nent between and 2. Notice that with this definition 
D* has units of distance square. 

Likewise, we can define the effective diffusion coeffi- 
cient for a subdiffusive and a superdiffusive limit regimes 
in our model. In the first case, setting e = and 7 = A 
in (|20|). we have that 



d: 



1 



ub 



s r(A + i) 



In the superdiffusive limit, setting 7 
(1^0)1 we obtain 



(27) 



2A — a in Eq. 



d: 



(2ry)2 



K 



X.a 1 



(28) 



where we have used relation (p4)l to derive this expres- 
sion. 

Interestingly, these analytical expressions enable a 
quantitative interpretation of the physics underlying the 
transport processes, linking the macroscopic effective dif- 
fusion coefhcient with microscopic parameters derived 
from the forces correlation functions. 

In Ref. [1^ we obtained the values of D* by fitting the 
MSD-vs-time lag plot of each trajectory with the empir- 
ical expression given by Eq. (j26p . The distributions of 
D* followed an exponential behavior, with mean values 
equal to 2 x 10~^/im^ and 4.4 x IQ'^jim^ for aggregation 
and dispersion conditions, respectively. 

Using ((28)) and the parameters displayed in Table 1 to 
compute the analytical effective diffusion coefficient, we 
obtained (5 ± 2) x lO'^^m^ and (4 ± 1) x lO^^^m^, for 
aggregation and dispersion, respectively. These values 
are consistent with the ones obtained empirically. 

For mutant cells, the computed analytical effective dif- 
fusion coeflicient ^ in D* ^ (0.85 ± 0.1) x 10-V"^^ 
which is in agreement with the one obtained in Ref. (Tsj 
(~ 1 x IQ-^^im^) 

Moreover, the correlation function ^ allows obtaining 
an estimation for the magnitude of the "global" force 
exerted by the motors Fmot, which is given by 



Frr 



An 



r(i-a) 



(29) 



An 



■(2r,)^ 



:{2kBTf 



(30) 



Taking hgT = 4 x lO-^pA^^m, {2'qf =MSDo and e, 5 
and a from Table I, we obtain Fmot 8 ± ApN and 
16±6piV, for dispersion and aggregation, respectively, 
which are in the order of the reported stall force of a 
single myosin-V j42i] . The values obtained for aggregation 
and dispersion are similar within the experimental error, 
since as far as we know the stall forces for myosin in these 
stimulation conditions are expected to be the same. 



VI. SUMMARY 

The violation of the FDT in living cells has been ob- 
served in recent studies [1, [ll|, [13, [lJ|- Motor proteins, 
which are force generators in cells, can not only mod- 
ify the viscoelastic response of the cytoplasm -as is the 
case of CSK motors-, but they also generate non-thermal 
random forces which drive the system out of equilibrium 
[m . The action of motors is generally reflected in a tran- 
sition between a subdiffusive to a superdiffusive behavior 
of the displacement of particles within the cells Q . 

To investigate how the action of molecular motors af- 
fects the transport in living cells, we study the motion 
of melanosomes driven by myosin-V along the F-actin 
network in Xenopus laevis melanocytes. 

Differently to other microrheology experiments re- 
ported before 0, H, [HI, [H, [i3l , we adopt a different ap- 
proach which uses a generalized Langevin equation that 
explicitly considers the collective action of the molecular 
motors and the presence of experimental/biological noise. 
The analytical solution of the model is obtained and an 
expression for the slope of the MSD vs the time lag is 
derived explicitly. The comparison between the model 
solution and the experimental data is straightforward, 
and involves 4 parameters. Two of these parameters are 
the exponents of the thermal and non-thermal forces cor- 
relation functions, while the others are related with the 
relative intensity between both kind of forces and the 
residual value of the MSD. 

On one hand, the model predicts the observed 
crossover between subdiffusive to superdiffusive regimes, 
as well as it gives good estimates for the in vivo motor 
forces. On the other hand, the proposed method en- 
ables us to determine a link between the macroscopic 
effective diffusion coefficient and the parameters in the 
microscopic scale. 

We believe that this theoretical approach can be used 
to describe the dynamics of intracellular transport of dif- 
ferent cargoes in other living cells. 



8 



Acknowledgments 



clonal de Promocion Cientifica y Tecnologica, Argentina. 



We acknowledge support from grants PICT 928/06, 
PICT 31980/05 and PICT 31975/05 from Agenda Na- 



R. Mallik and S.P. Gross, Curr. Biol. 14, R971 (2004); [23 
R. D. Vale, Cell 112, 467 (2003). [24 
S. P. Gross, M.A. Welte, S. M. Block and E. F. Wi- 
eschaus, J. Cell Biol. 148, 945 (2000). [25 
V. Levi, A. S. Serpinskaya, E. Gratton and V. I. Gelfand, 
Biophys. J. 90, 318 (2006). 

L. Bruno, M. M. Echarte and V. Levi, Cell Biochem. and [26 
Biophys. 52, 191 (2008). 

S. P. Gross, M.C. Tuma, S.W. Deacon, A.S. Serpinskaya, 
A.R. Reilein and V. I. Gelfand, J. Cell Biol. 156, 855 [27 
(2002). 

G. M. Langford, Traffic 3, 859 (2002); S. J. Atkinson, [28 
S. K. Doberstein and T.D. Pollard, Curr. Biol. 2, 326 
(1992). [29 
D. Mizuno, C. Tardin, C. F. Schmidt and F. C. MacK- 
intosh. Science 315, 370 (2007). 

P. Bursac, G. Lenormand, B. Fabry, M. Oliver, D. A. [30 
Weitz, V. Viasnoff, J. P. Butler and J. J. Fredberg, Nat. 
Mater. 4, 557 (2005). [31 
J. Snider, F. Lin, N. Zahedi, V. Rodionov, C.C. Yu, and 
S.P. Gross, Proc. Natl. Acad. Sci. USA 101, 13204 [32 
(2004). 

See T. A. Waigh, Rep. Prog. Phys. 68, 685 (2005) and [33 
reference therein. 

A. W. C. Lau, B. D. Hoffman, A. Davies, J. C. Crocker, [34 
and T. C. Lubensky, Phys. Rev. Lett 91, 198101 (2003). 
C. Wilhelm, Phys. Rev. Lett 101, 028101 (2008). 
M. Brunstein, L. Bruno, M. A. Despsito and V. Levi, [35 
submmited. 

F. Gallet, D. Arcizet, P. Bohcc and A. Richert, [36 
arXiv;0901.3087 (2009). 
C. Metzner, C. Raupach, D. Paranhos Zitterbart, and B. [37 
Fabry, Phys. Rev. E 76, 021925 (2007). 

L M. Kulic, A. E. X. Brown, H. Kim, C. Kural, B. Blehm, 

P. R. Selvin, P. C. Nelson, and V. I. Gelfand, Proc. Natl. [38 

Acad. Sci. USA 105, 10011 (2008). 

H. Salman, Y. Gil, R. Granek, and M. Elbaum, Chem. [39 
Phys. 284, 389 (2002). 

G. Lenormand, J. Chopin, P. Bursac, J. J. Fredberg, and [40 
J. P. Butler, Biochem. and Biophys. Res. Comm. 360, 

797 (2007). [41 
X. Trepat, G. Lenormand, and J. J. Fredberg, Soft Mat- 
ter 4, 1750 (2008). 

Y. Tseng, T. P. Kole, and D. Wirtz, Biophys. J. 83, 3162 [42 
(2002). 

M. Weiss, M. Eisner, F. Kartberg, and T. Nilssony, Bio- 
phys. J. 87, 3518 (2004). 

M. J. Saxton, Biophys. J. 86, 369 (2004); M. J. Saxton, 
Biophys. J. 92, 1178 (2007). 



D. S. Banks and C. Fradin, Biophys. J. 89, 2960 (2005). 
D. S. Martin, M. B. Forstner, and J. A. Kas, Biophys. J. 
83, 2109 (2002). 

G. T. Shubeita, S. L. Tran, J. Xu, M. Vershinin, S. Cer- 
meUi, S. L. Cotton, M. A. Welte, and S. P- Gross, Cell 
135, 1098 (2008) 

Y. Shtridelman, T. Cahyuti, B. Townsend, D. DeWitt 
and J. C. Macosko, Cell Biochem. Biophys. 52(1), 19 
(2008). 

A. A. Nascimento, J. T. Roland, and V. L Gelfand, Annu. 
Rev. Cell Dev. Biol. 19, 469 (2003). 

V. Levi, V. I. Gelfand, A. S. Serpinskaya, and E. Gratton, 
Biophys. J 90, L7 (2006). 

M. M. Rozdzial, and L. T. Haimo, Cell 47, 1061 (1986); 
P. J. Sammak, S. R. Adams, A. T. Harootunian, M. 
Schliwa, and R. Y. Tsien, J. Cell Biol. 117, 57 (1992). 
S. L. Rogers, L S. Tint, P.C. Fanapour, and V. L Gelfand, 
Proc. Natl. Acad. Sci. U S A 94, 3720 (1997). 
S. L. Rogers, R. L. Karcher, J. T. Roland, A. A. Minin, 
S. W., and V. L Gelfand. J. Cell Sci. 146, 1265 (1999). 
V. Levi, A. S. Serpinskaya, E. Gratton, and V. I. Gelfand, 
Biophys. J. 90, 318 (2006). 

C. Dietrich, B. Yang, T. Fujiwara, A. Kusumi, and K. 
Jacobson, Biophys. J. 82 , 274 (2002). 
P. Bursac, B. Fabry, X. Trepat, G. Lenormand, J. P. 
Butler, N. Wang, J. J. Fredberg and S. S. An, Biochem. 

and Biophys. Res. Comm. 355, 324 ( 2007). 

L. Bruno and M. A. Desposito, larXiv:0904.2T33V l 
[physics. bio-ph] . 

R. Zwanzig, Noneqmlibrium Statistical Mechanics (Ox- 
ford Univ. Press, New York, 2001). 

K. G. Wang, Phys. Rev. A 45, 833 (1992); E. Lutz, Eu- 
rophys. Lett. 54, 293 (2001); A. D. Vitiales and M. A. 
Desposito, Phys. Rev. E 73, 016111 (2006). 
S. Sharma, S. Wagh, and R. Govindarajan, Pigment Cell 
Res. 15, 127 (2002). 

S. Yamada, D. Wirtz, and S. C. Kuo, Biophys. J. 78, 
1736 (2000). 

P. Dieterich, R. Klages, R. Preuss and A. Schwab, Proc. 
Natl. Acad. Sci. USA, 105 , 459 (2008). 
C. Raupach, D. Paranhos Zitterbart, C. Mierke, C. Met- 
zner, F.A. MuUer, and B. Fabry, Phys. Rev. E 76, 011918 
(2007). 

G. Cappello, P. Pierobon, C. Symonds, L. Busoni, J. 
Christof, M. Gebhardt, M. Rief, and J. Prost, Proc. Natl. 
Acad. Sci. USA 104, 15328 (2007); A. E. M. Clemen, 
M. Vilfan, J. Jaud, J. Zhang, M. Barmann, and M. Rief, 
Biophys. J. 88, 4402 (2005). 



